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Abstract 



We discuss corrections to the linear response of a many-body system be- 
yond the binary cohision approximation. We first derive for smooth pair 
interactions an exact expression of the response oc \/q^, considerably simpli- 
fying existing forms and present also the generalization for interactions with 
a strong, short-range repulsion. We then apply the latter to the case of liquid 
^He. We display the numerical influence of the l/q^ correction around the 
quasi-elastic peak and in the low-intensity wings of the response, far from 
that peak. Finally we resolve an apparent contradiction in previous discus- 
sions around the fourth order cumulant expansion coefficient. Our results 
prove that the large-g response of liquid '^He can be accurately understood on 
the basis of a dynamical theory. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



High-energy pulsed neutrons from spallation sources have recently been used for the 
collection of good-quality cross sections data for the inclusive scattering of neutrons from 
liquid ^He. Data are for temperatures below and above the transition temperature Tc§'li. 
The above cross sections are a direct measure of the dynamic response or structure function 
S{q,uj), where g, a; are the momentum and energy transferred to the system. 

It appears convenient to consider the reduced response y) = {q/M)S{q, u) instead of 
S{q,uj), where the energy loss parameter u is replaced by an alternative kinematic variable 
y. The latter is a linear combination of {q,uj). M is the mass of a constituent atom. 

Virtually all dynamic calculations of the high-g response S(g,cu)Bi'ifll§ffl have been 
based on the Gersch- Rodriguez-Smith (GRS) expansion of the reduced response in 1/g or 
modifications of it@'iili'0. The theory in principle employs only the elementary atom-atom 
interaction V and is otherwise free of parameters. 

The dominant part of the large-g response is the asymptotic limit. It describes the 
response of a neutron striking an atom with given momentum. The absorption of the 
transferred momentum and energy-loss is not affected by other atoms in the medium. Final 
State Interactions (FSI) induced by V, produce corrections to the above limit which vanish 
for increasing q. The leading FSI oc 1/g is caused by binary collisions (BC) between a struck 
and an arbitrary second particle in the medium. For liquid ^He, predictions for the reduced 
response y) to order 1/g agree well with the data over a broad range around the central, 
quasi-elastic peak at y = 0. In fact, the quality of the data hardly calls for refinements 
beyond the BC. The incentive to nevertheless consider the introduction of fine details is 
mainly of theoretical nature: A criterion for the expansion of a function of two variables 
(f){q,y) in 1/g must depend on y. In particular the large \y\ wings where the response is 
only a small fraction of the peak value, has been suspected before to be sensitive to details 
beyond the BC. 

In this note we treat, to our knowledge for the first time, 1/g^ corrections. Those are 
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due to ternary collisions (TC) between a struck and two other particles. Their study is the 
major purpose of this note. 

A second topic to be discussed is related to the cumulant expansion of the response 
which has recently resulted in a successful, model-independent extraction of the single-atom 



apparent discrepancy between the directly computed fourth cumulant coefficient and the 
value extracted in the BC approximation for a dynamically calculated response0. 

We start in Section II recalling some essentials of the GRS expansion for the reduced 
response, valid for smooth inter-particle interaction! and derive a formally simple repre- 
sentation of TC terms oc Next we mention modifications which are required if the 
pair-interaction has a strong, short-range repulsion. In Section III we present numerical val- 
ues for TC contributions to the response of liquid ^He and discuss its relative importance, 
in the peak region and the low-intensity wings. Section IV contains a brief discussion of the 
cumulant representation of the response which has recently been used to parametrize data 
for liquid ^He. We report a complete calculation of the 4th order FSI cumulant coefficient, 
and thereby also resolve a previously existing discrepancy. 

II. DOMINANT FSI PARTS IN THE RESPONSE FOR SMOOTH AND 

SINGULAR INTERACTIONS. 

Consider the response S{q, uj) per particle for an infinite system of point-particles 



momentum distribution in liquid ^HeBByelBO. Our interest here is the resolution of 




an 



(2.)-'/ 



OO 



S{q,u)=A-^ 



rfte^-*($o|pJ(t)Pg(0)|$o) 



(1) 



— oo 



Above Pq{t) is the density operator, translated in time t by the Hamiltonian H 



iHt 



(2) 



$0 is the groundstate of the system with energy Eq. 
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We shall work with the reduced response (p{q,y) = {q/M)S{q,u!), where the energy loss 
00 is replaced by an alternative kinematic variable y = y{q,uj 

Substitution of (2) into (1) produces incoherent and coherent components. Considering 
high-g responses, it suffices to discuss the dominant incoherent part, where one tracks a 
single particle (for instance '1') in its propagation through the mediumEll. 

For the description of the large-g response we shall follow the formulation of Gersch, 
Rodriguez and Smith (GRS)S and cite from there a few results. Details can be found in the 
bibliography. 

It is convenient to introduce the Fourier Transform (FT) of the reduced (incoherent) 
response 

_ fOO 

(P{q,s)= dye-^^y<P{q,y) (4) 



The variable s = {q/M)t above is the distance traveled in the medium during a time t by a 
constituent, which moves with constant recoil velocity Vg = q/M: s is the length, canonically 
conjugate to the momentum y. 

The density fluctuations Pq in (Q) shift coordinates in the direction of q, chosen to be 
the z-direction. It leads to the following, standard expression, valid for local forcesSS 
(r — s = r — sq) 

4>{q, s) = (l/A)(^$o(ri - s; rfe)|T,exp|(Vwg) da[H{ri - a; - Eo]||$o(ri; r^)^ 

= E(lM)"^n(^) (5) 

n 

The second line in (5) is the GRS series in l/vq, which is generated by the expansion of the 
above, a-ordered exponential. 

For local interactions, the Hamiltonian with shifted coordinate '1 ' can be written as 

H{ri -cT-rk)=Y.Ti + Y. ^^'^h ru) + E ^('^i " ^k) 

I (7^1 fc>l 

k>l 



[Y^Ti+Y. V{rf, n)] + E [^(^1 - - V{r^; r,)] 

/ l,k>l k>l 

H{r,;n)+U,{a), (6) 



with 



f/i(a)(= f/i(a; n, r,)) = ^ 5.y(ri; r^) 

S.Vin; n) = [V{r^ - a; n) - V{r^; r,)] (7) 

SaV{ri; Vk) is the difference between the interaction of a selected particle k and '1' with the 
latter, once at a shifted position Vi — a and then at ri, Ui{a) is the same due to all particles 
A; 7^ 1. Using (P) one checks! 

[H{ri - a- r,,) - ^o]'^>o(?'i, r^) = Ui{a)%{r^, r^) (8) 

It is convenient to introduce the FT of the GRS coefficient functions in (^. For example 



-Fo(s) = ^ / dridrkPA{ri, ri - s; Vk) 



P 



^.-^P-^n{p) (9) 



(27r)3 

The dominant correction to the above asymptotic limit is (r = ri — 

1 ^ . . 2 



Vq ...q 



Fi{s) = J dridVkpAiru fi - s; r^) ^ d(r[H{ri - cr) - Eq] 
J dridrkpA{ri,ri - s;rk) daUi{a) 



AVq 

^./d /'<'--'-^^°> X.(r.4 (10) 

where use has been made of (^. Above the function x 

Xq{r,s) = -{l/vq) f da6,V{r), (11) 







is the off-shell, eikonal phase in the coordinate representation, pertinent to the characteristic 
difference of interactions acting on '1' in (0). In order to obtain the appropriate on-shell 
phase one needs to replace the integration limits in ( ]Tl| ) from 0, s to — oo, oo. 



In all, Fi{s) describes binary collisions (BC) of the struck particle '1' with any other 
particle in the medium. Its FT Fi{y) is odd in y and for large q mainly shifts the position 
of the maximum of the even asymptotic -Fo(z/) aX y = 0. 

The computation of the above quantities requires non-diagonal density matrices. We 
shall use a normalization, such that 

Pn{l, ...n; 1', ...n') = (^^f=n+i J ...n; n + 1, ..n^)$o(l', --n'; n + 1, ...ua) 

= [{A-n- 1)]-^ J d[n + l]p„,+i(l, ...n, n + 1; 1', ...n', n + 1) 

Pa{1, ...A; 1', ...A') = A\%{1, ...A)%{1', ...A'), (12) 

The densities p„(ri,ri — s^r^), required in Eqs. (|), (p^Of ) , are diagonal in all particles 
k, except in '1'. For example, pi(0;s) = pi{r,r — s) in (|) is the single-particle density 
matrix which has the single-particle momentum distribution n{p) as its FT. Its diagonal 
part pi(0, 0) = p is the number density. 

The first two terms in the GRS series Fq, Fi satisfactorily describe the data for the 
reduced response in a broad band around the quasi-elastic peak at y = 0. Such a fit 
cannot be expected in the wings, where Fi{y) occasionally reaches small negative values. 
For growing y in those wings Fi{y), I > 2 competes with F^^y) + {l/vq)Fi{y) of comparable 
size: ultimately 4>{q,y) > 0. 

Partial expressions for the next-to-leading order terms F2{s), F2{y) have been given before 
by Gersch et ali and by Besprosvanylli but those forms are not complete and are not always 
transparent. We shall derive below expressions for the exact F2, based on Eqs. (12), (15) in 
Ref. I 

his) = -jJ dridrkpA{ri, ri - s; rt) da[H{ri - cr, Vk) - ^0] ^ da'[H{ri - cr', r^) - ^0] (13) 

Consider the operators in the brackets above, acting on the wave functions which compose 
Pa, Eq. (p^, with shifted, respectively unshifted coordinate ri. Their combined result is 

F2is) = -^J dridVkPAiri, Vi - s; r^) da\Ui{a) - Ui{s)] da'Ui{a') 

= jJdr,drkpA{ruri-s;rk){^ daUi{a) ~ Ui{s) j'^ da da'Ui{<y')^ (14) 



Since F2 oc f/^ and the latter oc V'^ we dub F2 the TC contribution to FSO 

We start with the first term in the braces in Eq. (p!^. Using the definition in one 
rewrites 



J2 [ da6^V{ri - 
+ r da6„V{ri-rk) da 6^V{ri - Vi) 



(15) 



The above components are distinct two- and three-particle operators and the same holds 
for the parallel decomposition of the second term in the braces in (|14D- With ri = Vi — 
^2; ^'s = ■'"3 — 1^2) one may reduce ([T^) to 



V 



F2M 



dr 



P2{r - s,0;r,0) 



P 



\[iXq{r, s)f + -6sV{r) f da[iXq{r, a)] 

I Vq Jo 



(16a) 



dvidv; 



Pzjri - 5,0,r3;ri,0,r3) 
P 

+—6sV{ri) da[ix{ri-f3,a)] 
v„ Jo 



■MXgiri,s)][iXqiri -r3,s)] 



(16b) 



with Xg(^) s) the off-shell phase as defined in (|TT]) (Note that x 1/^g)- 

For smooth, non-singular local forces, the above completes the derivation of an exact 
expression for F2{s). However, if V possesses a strong, short-range repulsion, as is the case 
for atom-atom interactions, difficulties emerge. There are no problems if in integrands wave 
functions or density matrices and V have identical arguments, in which case large values of 
V are generally off-set by small values of p„ at common small r. This is not the case in 
general. A prime example is Eq. ([I0|) with r-dependence through V{r — cr), < |<t| < s, 
and r — <T generally not coinciding with either r — s or r in p2{r — s,r;0): Large line 
integrals may result. 

In the above case smooth expressions emerge again upon summation of a ladder of pair 
interactions V{r), leading to VJ=// which is the eikonal, off-shell t matrixi'i0. Effectively 



iX exp[zx] - l = ix + l/2[^x] + •• 



(17) 



Using ([1^) and ([TTD we define 



G2{s, [V]) ^ F2{s,V ^ [t]) 



(18) 



with the following two- and three-particle components 



V 



1 



dr 



P2{r - s,0;r,0) 



P 



—SsV{r) da[exp[ixg{r,a)]-l 



(19a) 



P3{ri - s,0,r3;ri,0,r3) 
dridr^ 

P 



exp[zxy(ri - rg, s)] - l) + ^SsV{ri) da(^exp[ixg{ri - rg, a)] - 1 



(19b) 



Care should be exercised in the replacement V — Ve// for an ill-behaved V as we shall now 
illustrate by focusing on the first term in the brackets in ( |16a| ), It can be shown that 

all higher order terms Fn{s, [V]) contain a 2-body component of the form 



^Fi'\s, [V]) 



dr 



P2{r - s,0;r,0) 



{ixr/n\ 



(20) 



Using (IT3) those may be summed up to 



(21) 



jt>i 



If therefore Fi{s, \V]) has been regularized by — K// = t, the first part of F2 (s) in Eq. 



( |16a|) is already contained in -Fi(s, [t]) and in order to avoid double-counting it should be 
removed from 6*2 (s, \y])- We note that by construction, the remaining term in Eq. (|19a|) 
which emerges from TC is nevertheless of 2-body character. 

Eqs. (|19D above are the exact TC contributions to the GRS series (|^) for the reduced 
response. At this point we mention an alternative to the GRS series, namely the cumulant 
representation for the FT of the reduced responsel^ 



0(g, s) = Fo{s)R{q, s) = Fo(s)exp[n(g, s)] 



n>l ^'^1 



Fn(s) 



h={R-l)-il/2)iR-lY + ..., 



(22a) 
(22b) 
(22c) 
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with all FSI effects contained in either i? or f2 

Below we report a calculation of TC contributions, using choices for the underlying 
densities. 



III. TC CONTRIBUTIONS TO THE RESPONSE OF LIQUID "^HE. 

Until now, dynamical calculations based on the GRS series were limited to FSI in the 
BC approximation, i.e. 

n n^^ = n2 = - 1 (23) 

VqFo 



with the GRS series, cut at = 2, as in Ref. 0. We now report what apparently are 
the first results for the next-to-leading order TC corrections and which are contained in 
G2 = Gf + Gf ,Eqs. ®, dH). 

We first recall the standard input described in Ref. |rT[ For the bare interaction we use 
the standard V^"^, Ref. |2l], and for the single-atom momentum distribution the results 



of Refs. As regards semi-diagonal 2-particle density matrix p2{r — s,0;r,0), there 

exist results obtained using stochastic methodsBi, but computationally it is unnecessarily 
time-consuming to evaluate those for each and every (r, s), as required in calculations of the 
expressions (11^) or ( p.9a| ). 

In the past relatively simple guesses have been made for p2- We shall use below the 
interpolation formula by GRsi'i 

P2{r - s, 0; r, 0)(= p2{r - s, r; 0)) = ppi(0, s)C2{r - s, r) 

C2{r - s, r) ^ ^ g{\r - s\)g{r)) (24) 

with g{r) the pair-distribution function, chosen to be the one from Ref. ^. 

A calculation of G3 requires the 3-particle density matrix pa which, as before is non- 
diagonal in coordinate 1. As an approximation we suggest 

n - - n-\ (^-2)p2(ri-s,0;ri,0)p2(r3-s,0;ri-f3,0) 

P3 ri - s, 0, ra; ri, 0, r^) ^ — — 

[A-l) Pi(0;s) 
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= ^^P^P^(^'^ ^)C2(^i - ^•i)C2(ri - rg - s, f-i - ra) (25) 

where use has been made of (0). The choice ( pS]) has several advantages 

i) With 'V paying a special role, it is symmetric in the other coordinates 

ii) It exactly satisfies the 'sumrule' ( |I2| ) 

iii) It factorizes in parts dependent on fi,ri — 

iv) It causes pa to vanish for small values of the 4 coordinates which would otherwise 
produce large values for the factors in the operator in the brackets in (|19b| ). 

An immediate consequence of iii) above is the reduction of the, effectively 5-dimensional 
integral in ( |19b|) to the product of two, 2-dimensional integrals 



1 Gf\s) 

1 Gf{s) 
^1 h{s) 



P I dr(2[r - s,r 
-1 



P 



J driC2{ri - s, fi) (^exp[ixq{fi, s)] - 1 



(26a) 

J driC,2{ri - s,ri)5sV{ri) 

(26b) 



Anticipating small TC corrections we approximate the cumulant representation (|2^) 



R 



TC _ i^TC 
G2 



TC 



(27) 



The thus defined TC contribution to the FSI phase has been added to the previously calcu- 
lated BC part Cl^'-^ , Eq. (pS]). From ( P2a| ) and the inverse of (H), we compute the response 
for T = 2.5 K to the corresponding order. 

A first observation is the relative insignificance of 3-body TC contributions for the g-range 
investigated. A heuristic argument runs as follows. If the BC FSI contributions amounts to 
a fraction of the the dominant asymptotic limit, one estimates from the factorization of 
the 3-body density matrix that TC FSI is approximately the square of that fraction of Fq. 

We now display some results for TC contributions. Figs. la,b show for small y and 
q = 21, 25, 29, 50, 100 A ^ the even part of the calculated reduced response 0^^'^"(g, ?/) = 
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.y) + (j){q, —y)]/2 , without and including TC contributions (note that 0^*" is even). 
Even for y = there is an effect which for increasing g > 2lA ^ decreases from 2to0%. 
Figs. 2a,b show the fractional effect of TC contributions 

calculated for 2.5 ^ \y \ (in A ^) ^ 3.3. The difference in sign of a — 1 clearly shows the effect 
of the competition between the dominant even and odd parts in the wings of the response. 

Finally, Figs. 3a,b,c show for q = 21,25,29A~\t = 2.5K the efffect of TC on the 
calulated response, including the unresolved effect of the instrumental resolution. The small 
TC contributions discernibly improve the agreement of predictions with data in the above 
y regions. 

IV. ON THE FOURTH CENTRAL MOMENT OF THE RESPONSE. 



The preceding Sections deal with the reduced response ( [2^ ) up to, and including TC con- 
tributions and its, in principle exact, calculation. For the FT of those one needs F„(s), n < 2 
or (cf. ^2])) Q{q, s) to that order, both for all relevant s. 

We now address a second topic which is related to the cumulant representations (^21) and 
which is based on the small-s expansions 

Ms) = E —^'^m (29a) 
ms)=E-^M (29b) 
The above coefficients am are related to even moments of the momentum distribution n{p), 



while the FSI coefficient functions fSmiq) in the expansion (|29b| ) can be expressed in terms 



of central moments of the response (see for instance Ref. ^5]) 



Mniq)= j duo{uo-q^/2MYS{q,uo) 

= {v,r I dyy-ct^{q,y) ^ (t;,)"A<„(g) (30) 
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For our purpose it suffices to give the following expressions for n=3,4 and valid for local 
interactions ^0 

^3 = -M3= (^)(V^V^) (31a) 
P^ = M^-a,- Ml = {J'l.J^i) (31b) 



J-\ above 

J=r = Y. ^) = E ^(^1 - ^A.), (32) 
fc^l fc>l 

distinct from ?7i(s), Eq. (|^, is the true total force on a given particle 'I'll. The expectation 

value in (|31b| ) can thus be separated in two parts. The first contains the square of the force 

on '1' due to one particle and in the second part forces on '1' by two different particles 

(^i-^i) = (E[^i(l>j)]') + ( E :^i(l,j)^i(l,^)) (33) 



The expansions (|29|) provide a parametrization of the response, but the technique has been 
shown to have its problems0. One such problem is the convergence for growing s which is 
indispensable for the calculation of the inverse FT (^) from (j){q,y). Moreover, the cumulant 
expansion lacks a systematic ordering in powers of 1/g which is also remedied in GRS 
theory. Notwithstanding, there has recently been a renewed interest in the above small 
s-cumulant expansions as a vehicle to extract the single-atom momentum distribution n{p) 
from response data for ^He and NefflHilffl. Around an apparent contradiction arises, 
which we discuss below. 

One may 'invert' Eqs. (|2^) in order to find alternative expressions for 



^^(g)=m!^™ \imKl{q,s)/s' 



In particular 



/33(g) = 6 lim [Imil(g, s)/s^] (34a) 
^4(g) = 24 lim [Refi(g, s)//] (34b) 

s— >0 
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The above FSI coefficient functions can only be calculated if a theory provides the FSI phase 
function Cl{q, s). The GRS theory is one such example. The dynamic calculation described 
in the previous Sections, provides Q{q, s) for all s. 

First we state that without truncations, the cumulant expansion and the GRS series 
ought to lead to the same response, and in particular to the same numerical values for the 
cumulant coefficient functions Pm- This is not self-evident since Eqs. (|3l|) and (^) look 
quite dissimilar. The former are expectation values in terms of diagonal density matrices, 
whereas the Q underlying the GRS theory is an operator, averaged over a non-diagonal 
density matrix. Nevertheless the identity of the derived Pslq) has been formally verified in 
the past (see for instance Ref. |^). In addition a numerical test has been performed using 
Q — s> Cl^'-' , which suffices since (3^ draws entirely on Fi or on Cl^*-^ . The calculated value and 
the one, extracted over a wide q range, indeed agree to high accuracy0. 

For a similar demonstration regarding ^4 one uses (|22b| ), (0), ([T5| ) in ( P4a| ) and readily 
verifies that terms oc s^, needed in the threshold behaviour ( |34b| ), originate exclusively from 
the TC terms G2{s) (cf. Eqs. (0)). Observing that 6V in (13), (14) always appears 
quadratically, one has 

UU (X s\z.:fI + 0{s^), (35) 

and ( p^ ) results. 

We separately treat BC and TC contributions to ^4(13') and start with the above men- 
tioned BC approximation for the regularized Fi{y, [t]). One observes that even the BC FSI 
phase function Cl^'-^{q,s) contains terms oc s^, contributing to f]^. Again, a remarkably 
stable, negative value could be extracted from calculated BC phases over a wide g-range0 
{q* = q/10 in A~^) 

g*2^f 1) (g) = (-2.27 ± 0.02) A-^ (36) 

One easily demonstrates that the same, computed from the ffist term in the brackets of 
( |l6a|) is g*2;5f^)(g) = -(M/IO)^ S^.^i(:Fi(1, j)2) = -2.19A-^ again in close agreement 
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with the extracted result (pGl). The negative outcome clearly contradicts the manifestly 
positive expression ( p 1 b| ) for the complete P^. The latter, however, draws also on additional 

TC contributions from and G^2 \ Eqs. (p!9|), which we now address. 

~ (2) 

We start with the threshold value of the two-body part G2 of the TC contribution, 
which is readily shown to be exactly 4/3 times the first part and for the positive, complete 
two-body part one finds (cf. (^3|) ) 



(g) = q*'[Pt\q) + Pt\q)] = l:iM/ioYJ2{^^ii,jf) = 0.73A-' (37) 
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Within ^ 0.5% the same value results when calculating the threshold value ( |36| ) and using 
for the FSI phase function il{q,s) Eq. (p2D with as in ([l8|) , ( |19a| ). As emphasized 



before, the close agreement evidences numerical accuracy and not consistency. 
The genuine 3-body TC part, defined by (|33D , |19b|) 



^i^\q) = 241im ( — ) Re 



*-»0 \Vq 



1 



2 



Fnis] 



(38a) 



= ^(-)'(M/10)^ E (^i(l,j)-:^i(l,A:)) (38b) 
involves the forces on '1' by two different medium particles. The expectation value in 



requires a diagonal 3-particle density matrix, and consistency requires it to be the non- 
diagonal (^), used in the calculation of Q in the limit s = 0. Here too the 3-body part is 
neghgible. 

The remaining 2-body parts may be compared with previous results which have been 
calculated in different ways. Our result (p6D lies in between (0.69, 0.86) A"^, communicated 
by Polls from various approximations to the pair-distribution functions Another 



stochastic calculation by Clyde and Boninsegni, reported in Ref. |T4|, leads to a result about 
35% in excess of the above. Previous experience has taught that averages, like the ones in 
( piap and ( |31b| ) are quite sensitive to the chosen, pair distribution. A spread of 10-15% may 
certainly be expected, but presumably not a deviation of ~ 35%. 

Finally, we compare the computed total q*^Pi{q) = 1.17A~^ with a few results, extracted 



from cumulant analyses of the data. For instance in Ref. ^ a value compatible with is 
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given, while an upper limit q*'^Pi{q) < O.SOA"'* is cited in Ref. |6[ 

We close this Section by comparing F2{s), Eq. (p!^, and other published expressions for 
the sameBislS. Those are also quadratic in V , but contain in addition to p2, derivatives 
of p2 and V . In contra-distinction our result is quadratic in V and free of derivatives. We 
have shown above that each of its composing parts is oc [^^] with different co-factors. The 
alternative expressions provide directly one factor !Fi and it is not at all evident that the 
other part can be cast in that form. The equations of motion for density matrices ultimately 
provide the evidence. The procedure followed in Eqs. (3)-(6) avoids those steps and leads 
directly to the desired result. This can be checked for the general response of a particle in 
a potential, Eq. (8c) of Ref. ^ 



V. SUMMARY AND CONCLUSION. 

We have derived above an exact expression for the contribution of ternary collisions to 
the response of a non-relativistic many-body system, where the struck constituent interacts 
with two other medium particles. Its numerical contribution has for the first time been 
evaluated for the response of liquid ^He, T > and momentum transfers in excess of 
21 A ^. For those we know that the asymptotic limit and the dominant binary collision 
correction, accurately describe the response in a broad region around the quasi-elastic peak, 
but not necessarily at the peak itself (cf. for instance Ref. [11]). 

Our main interest was therefore focused on y ^ and the region of the wings, where the 
intensity is only a fraction of that in the peak. Compromising only on the assumed 3-body 
density matrix, we computed the relative size of small TC FSI effects and found those to 
discernibly improve the agreement with the data. 

The above calculation completes a program to calculate the medium-to-large q response 
of liquid ^He. A number of conclusions are in order. Using exclusively the well-known 
atom-atom interaction, basic ground-state properties as are the single-atom momentum 
distribution, the pair-correlation function and non-diagonal, two-particle density distribution 
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have been determined with great accuracy. 

The above quantities are then basic input for the calculation of the linear response of 
the system. Only weak assumptions have been used for the required two- and three-particle 
density matrices, diagonal in all, except one coordinate. Excellent agreement has been 
obtained with data for a theory with demonstrated convergence. 

Indeed, given the non-negligible scatter in the data and observing that one deals with 
atomic dynamics and not with QED, we feel that there is at present no incentive to study 
even finer theoretical details than discussed up to now. 

Our final remark regards the response of liquid "^He when compared with the responses 
of other systems, composed of atoms, molecules, atomic nuclei or sub-hadronic matter. We 
do not know of a system where the approach to the asymptotic limit has been measured 
and studied with an accuracy, possible for ''He. 
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Figure Captions 

Figs. la,b. The approach to the asymptotic hmit FQ{y) (diamonds) for small y of the 
calculated even part of the response without and with the (even) TC contributions. 

Figs. 2a,b The fractional effect a{q,y), Eq. (p^) of TC contributions in the wings 
2.5 < |i/|(inA"') <3.3. 

Fig. 3a. Calculated response and data for q = 2lA ^ , T = 2.5 K, including the effect of 
instrumental resolution. Dashed and drawn curves are without, respectively including TC 
contributions. 

Fig. 3b. Same as Fig. 3a for q = 25 A ^ . 

Fig. 3c. Same as Fig. 3a for q = 29 A ^ . 
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